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Abstract. We present 2D axisymmetric and 3D numerical simulations of the proto-planetary nebula Hen 3-1475, which is 
characterized by a remarkably highly coUimated optical jet, formed by a string of shock-excited knots along the axis of the 
nebula. It has recently been suggested that the kinematical and morphological properties of the Hen 3-1475 jet could be the 
result of an ejection variability of the central source (Riera et al. 2003). The observations suggest a periodic variability of 
the ejection velocity superimposed on a smoothly increasing ejection velocity ramp. From our numerical simulations, we 
have obtained intensity maps (for different optical emission lines) and position-velocity diagrams, in order to make a direct 
comparison with the HST observations of this object. Our numerical study allows us to conclude that a model of a precessing 
jet with a time-dependent ejection velocity, which is propagating into an ISM previously perturbed by an AGB wind, can 
succesfully explain both the morphological and the kinematical characteristics of this proto-planetary nebula. 
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1. Introduction 

Recent images of young PNe obtained with the HST have re- 
vealed that almost all of them have aspherical morphologies. 
Also, a significant fraction show bipolar or mutipolar lobes, 
highly collimated outflows or jets and other complex structures 
such as loops (see, e.g., Sahai & Trauger 1998; Garcia-Lario, 
Riera & Manchado 1999; Trammell & Goodrich 2002). Sahai 
& Trauger (1998) suggested that the outflows (i.e., collimated 
fast winds, hereafter CFWs) acting during the early post- AGB 
phase may have an important impact in the formation of these 
complex structures. However, the mechanism(s) producing the 
collimated fast winds or jets are not well understood. 

The most promising scenarios for producing CFWs or jets 
are : 

(i) magnetized accretion disk and binary stellar system (see, 
e.g., Soker & Livio 1994; Blackman et al. 2001), 

(ii) a rapidly rotating star with a toroidal magnetic field (see, 
e.g., Garcia-Segura et al. 1999; Garcia-Segura & Lopez 2000). 
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* Based on observations made with the Hubble Space Telescope, 
obtained from the Data Archive at the Space Telescope Science 
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Recently, the interaction of CFWs and the AGB envelope 
has been explored both analytically and numerically (Soker 
& Rappaport 2000; Lee & Sahai 2003; Garcfa-Arredondo & 
Frank 2003). 

The proto-planetary nebula (PPN) Hen 3-1475 shows a 
very well collimated jet, formed by a string of knots distributed 
along the main axis of the nebula (reminiscent of a HH ob- 
ject). Hen 3-1475 has been the object of optical imaging from 
the ground (Riera et al. 1995) and with the Hubble Space 
Telescope (HST) (Bobrowsky et al. 1995; Borkowski, Blondin 
& Harrington 1997; Ueta, Meixner & Bobrowsky 2000), near- 
infrared photometry (Garcia-Lario et al. 1997; Rodrigues et 
al. 2003), optical polarization (Rodrigues et al. 2003), opti- 
cal spectroscopy (Riera et al. 1995; Bobrowsky et al. 1995; 
Borkowski & Harrington 2001; Sanchez Contreras & Sahai 
2001; Riera et al. 2002, 2003), and proper motion studies 
(Borkowski & Harrington 2001; Riera et al. 2002, 2003). 

The central region is probably a dusty disk (or torus). The 
densest parts of the disk are seen in the WFPC2 images of 
Borkowski et al. (1997) as dark patches. The intrinsic polar- 
ization and large reddening of the circumstellar region of Hen 
3-1475 has been studied in detail by Rodrigues et al. (2003). 
The presence of CO molecular emission suggests the existence 
of large amounts of neutral gas. The CO emission shows an 
expanding torus of material (Knapp et al. 1995, Bujarrabal et 
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al. 2001). The recent CO observations of Huggins et al. (2004) 
show the presence of an expanding biconical envelope, and pro- 
vide strong evidence for entrainment of the molecular gas by 
the high-velocity jets. The OH emission suggests the existence 
of an aspherical shell expanding at 25 km s ' (Bobrowsky et 
al. 1995, Zijlstra et al. 2001). 

The radio continuum indicates the presence of a compact, 
ionized region surrounding the central B star (Bobrowsky et al. 
1995; Knapp et al. 1995). Optical, ground based spectroscopy 
of the central source revelead the presence of P-Cygni Balmer 
line profiles (Riera et al. 1995), which indicates the presence of 
mass outflows at high velocities. An analysis of a STIS, long- 
slit Ho- spectrum of the central region of Hen 3-1475 (Sanchez 
Contreras & Sahai 2001), suggests the presence of two different 
winds outflowing from the central star at high velocities. The 
highest velocity outflowing wind is highly coUimated close to 
the central star, and shows a radially increasing velocity. 

The emitting knots show extremely broad, double-peaked 
emission line profiles (with FWHM from 400 to 1000 km s ', 
Borkowski & Harrington 2001; Riera et al. 2003). The radial 
velocities observed along the Hen 3-1475 jet decrease with in- 
creasing distance to the central source, showing abrupt velocity 
changes in the transition regions from each (sub)condensation 
to the next one. 

From a comparison of the spectra observed for the interme- 
diate knots of Hen 3-1475 with the spectra predicted by plane- 
parallel shocks, Riera et al. (1995) concluded that the emis- 
sion is produced in a shock wave which propagates through 
a nitrogen-enriched medium. Predictions from plane-parallel 
shock models (adopting Type I FN abundances) with shock ve- 
locities from 100 —> 150 km s"' qualitatively reproduce the 
observed spectra (Riera et al. 2003). 

Recently, Chandra has detected X-ray emission from Hen 
3-1475 (Guerrero et al. 2004; Sahai et al. 2004). The X- 
ray source position coincides with the innermost optical knot 
NWl, and the X-ray emitting gas is probably shock heated. 
The X-ray spectrum implies a temperature of ~ 2.4 x 10'' K, 
and the observed X-ray luminosity implies a shock velocity ~ 
400 km s-i (Guerrero et al. 2004; Sahai et al. 2004). 

Recently, it has been suggested that kinematic and morpho- 
logical properties of the Hen 3-1475 jet could be the result of 
a variability of the central source (Riera et al. 2003). During 
the past 15 years, time-dependent ejection velocity models 
have been studied and applied to HH jets, and these mod- 
els in principle can also be applied to the Hen 3-1475 jet. In 
these models, the chain of knots observed in the jet is assumed 
to correspond to the successive internal working surfaces' re- 
sulting from a (possibly periodic) ejection velocity variability. 
One-dimensional analytic (Raga et al. 1990; Canto, Raga & 
D'Alessio 2000), as well as axisymmetric and 3-D numerical 
simulations (VoUcer et al. 1999; Raga et al. 2002; Masciadri et 

' A "working surface" is a double shock structure produced by the 
propagation of a fast flow into a slower flow. This working surface is 
formed by a bow shock sweeping up the material of the slow flow, an 
inward shock or Mach disk that decelerates the fast flow gas, and a 
contact discontinuity between them, separating the shocked slow and 
fast flow gas 



al. 2002) have been carried out for different functional forms 
for the ejection velocity variability. 

In the present paper, we investigate whether or not the prop- 
erties of the jets of Hen 3-1475 can be the result of the variabil- 
ity of the ejection velocity. As far as we are aware, our work 
is the first attempt to obtain a gasdynamical simulation of this 
object. The origin of both the collimation mechanism and of 
the ejection velocity variability is not studied in this paper 

The present paper is organized as follows : in section 2, a 
variable velocity ejection law is obtained for the case of Hen 
3-1475; the initial conditions and assumptions for the 2-D and 
3-D numerical simulations are presented in section 3; in sec- 
tion 4, we present the results of the temporal evolution from 
2-D axisymmetric numerical simulation (subsection 4.1) and 
the comparison of 3-D results with HST observations (subsec- 
tion 4.2); finally, in section 5, we give our conclusions. 

2. A variable velocity ejection law for Hen 3-1475 

We now consider the ground-based and the STIS HST spec- 
troscopy, together with proper motion measurements (taken 
from Riera et al. 2003) of the Hen 3-1475 jet. We first depro- 
ject the distance z from the source and the velocity v{z) for an 
assumed angle <p - 50° between the outflow axis and the plane 
of the sky (Borkowski & Harrington 2001). In this way, we ob- 
tain the physical velocity of the emitting gas as a function of 
distance from the source (see Figure[l). 

From Figure^ we see that the flow velocity decreases in an 
approximately monotonic way as a function of increasing dis- 
tances from the source. Such a velocity configuration suggests 
a time-dependent ejection velocity which smoothly increases as 
a function of time (with larger velocities at more recent times). 
Also, we need to superpose (on this "smooth ramp" variability) 
a variation in the ejection velocity with an appropriate period 
and amplitude in order to obtain internal working surfaces at 
positions that approximately correspond to the positions of the 
knots along the observed jet. 

After trying a number of different functional forms for the 
ejection velocity variability, we have chosen the form : 

Vj — vo + vi sin (2 7r(f - fo) / t) + at, (1) 

which includes a periodic term (described by a sinusoidal func- 
tion characterized by an amplitude velocity vi and a period r), 
superimposed on a linear increase of the velocity (at a rate a). 
In Eq. (1), t()(< 0) is the time at which the source is initially 
turned on (f = corresponding to the time at which the ob- 
servations were obtained). Through the computation of several 
models with different parameters, we find that the following pa- 
rameter set : V() - 400 km s"' (mean velocity), vi= 150 km s"' 
(velocity variability half-amplitude), t - 120 yr (period of the 
sinusoidal veriability), a - 1 km s ' yr ' (slope of linear ramp) 
and to - -640 yr (time at which the jet is turned on) produces 
jets with knot structures and kinematical properties similar to 
the ones of the Hen 3-1475 jet. We then use the ejection ve- 
locity variability given by equation (1) with this particular set 
of parameters in order to compute the axisymmetric and 3D jet 
simulations discussed in the following sections. 
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Fig. 1. Velocity of the emiting gas as a function of distance z 
from the source, reconstruted from the ground-based radial ve- 
locities (solid circles) and the STIS spectroscopy (open circles) 

3. Numerical Simulations 

3.1. 2-D Numerical simulation 

The 2-D numerical simulations have been carried out with 
the axisymmetric version of the Yguazii-a code (which is de- 
scribed in detail by Raga et al. 2000). This code integrates the 
cylindrical gasdynamic equations (employing the "flux vec- 
tor splitting" scheme of van Leer 1982) together with a sys- 
tem of rate equations for atomic/ionic species. With these 
rate equations, a non-equilibrium cooling function is com- 
puted. The reaction and cooling rates are described in detail 
by Raga et. al (2002)| We have used the mean Type I PN abun- 
dances from Kingsburgh & Barlow (1994). 

The axisymmetric numerical simulations were computed 
on a 4-level binary adaptive grid, with a maximum resolution 
of 1.172x10^^ cm along both axes. The computational domain 
is of (12.0, 1.5)x 10'^ cm, along the axial (z) and radial (r) 
axes, respectively. 

The jet is injected at z = with the time-dependent ejec- 
tion velocity given by Eq.([l) (and with the paremeter set given 
in §2). The jet has an initial radius ij - 6 x 10'^ cm, (time- 
independent) density of 500 cm"^ (similar to the density de- 
rived by Riera et al. in preparation) and temperature T^ = 1000 
K. 

We have carried out 2-D numerical simulations with two 
density distributions for the surrounding medium. First, we 
have assumed that the surrounding environment is homoge- 
neous, with a density n;.,,,, = 10 cm""* and temperature T;.,,,, = 
100 K. 

However, the high-velocity outflows of Hen 3-1475 are 
moving through the slow wind previously ejected by the star 
as a red giant on the AGB. The circumstellar envelope of Hen 
3-1475 is far from being spherical as shown from the optical 
polarization maps and the spatial distribution of the CO and 
OH emission (see, e.g., Bobrowsky et al. 1995, Zijlstra et al. 
2001, Bujarrabal et al. 2001, Rodrigues et al. 2003, Huggins 
et al. 2004). In order to include this, we have computed a sec- 
ond model, in which the outflow is moving into an AGB wind. 
For this model, we have assumed that the AGB remnant has an 
aspherical density distribution with high densities on the equa- 
tor and lower densities on the poles, which is described by the 
following equation (taken from Mellema 1995), 



P(R) = Po 



1 - a 



1 - expi—lfisirP'ff) 
1 - expi-ip) 



(Ro/Rf 



(2) 



where R - Vr^ + and 6 is the angle with respect to r-axis. 
The parameter a determines the ratio between the density at the 
equator and at the pole, and /3 determines the way the density 
varies from the equator to the pole (see Mellema 1995). The 
value of Po is calculated from the mass loss rate as : pQ - M 
I A n Venv, where Vf,,,, is the expansion velocity of the AGB 
remnant. We have chosen i?o = r^, a mass loss rate of 10 M0 



Fig. 2. Time evolution of the density stratification for the 
case of an outflow moving into a homogeneous environment. 
Integration times range from 100 to 600 years (top to bot- 
tom panels, respectively). The logarithmic greyscale, in units 
of g cm"^, is given by the vertical bar on the right of the bot- 
tom panel. The scale in both axes is given in units of 10'^ cm. 

Fig. 3. Same as Fig|2]but for the case of a jet expanding into 
an AGB remnant with a density distribution given by Eq|2](see 
text) 

yr"' and an expansion velocity of 20 km s"' . For a and (3 we 
adopt values of 0.7 and 3.0, respectively. With these parameters 
and Eq.(|3 the range of the AGB wind density along z-azix is 
~ [500, 0.05] cm"^ from the injection point to the z boundary, 
respectively. 

3.2. 3-D Numerical simulation 

Axisymmetric numerical simulations are powerful tools for 
finding a time-dependent ejection velocity (such as the one 
given by Eq.Q) which produces knots with appropriate po- 
sitions and kinematical properties. However, the Hen 3-1475 
jet has a "S"-like morphology, which is far from axisymmetric. 
Then, in order to reproduce this "S"-like shape, we have car- 
ried out 3-D numerical simulations, including a precession of 
the outflow axis with a period of 1500 yr and a half aperture 
angle for the precession cone of 7.5°. This precession produces 
approximately one half "wave" of the approximately helical jet 
beam within the observed length of the jet in Hen 3-1475 (in 
agreement with the observed images). 

We then carry out a 3-D numerical simulation of a jet with 
the ejection velocity variability given in Section 2, and the 
precession described above. The numerical simulations were 
carried out with the 3-D version of Yguazu-a code, employ- 
ing a five-level binary adaptive grid with a maximum resolu- 
tion of 2.34x10''' cm. The computational domain was of (3, 3, 
12)xl0'^ cm along x-, y- and z- axes, respectively. 

The outflow is injected at jc = y = 1.5 x 10'^ cm, on the 
z = plane. The initial jet radius is 1.5x10'* cm, the initial 
number density is ny = 500 cm"^ and the temperature is Tj 
- 1000 K. We have adopted a larger initial jet radius than in 
the axisymmetric simulations (see §3.1) in order to resolve the 
jet beam with ~ 10 grid points (in the lower resolution, 3-D 
simulation). 

In the 3-D simulation, we assume that the jet is moving 
into the AGB wind, for which we adopt the aspherical density 
function and parameters described in Section 3.1. 

4. Results 

4.1. Temporal evolution of jet structure for 2D 
numerical simulations 

Figure|2]shows a 600 yr time-sequence of the density stratifica- 
tion (in the rz-plane) for the case of a jet moving into a homo- 
geneous environment. Figure^shows a similar time-sequence. 
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Fig. 4. Comparison between simulated Ha emission for a 
model with a jet propagating into an uniform ISM (bottom 
panel), and the corresponding emission for the case of a jet 
moving into an AGB wind (top panel). Both maps have been 
obtained integrating the Ha emission coefficient along lines of 
sight, for an angle of 50° between the flow axis and the plane of 
the sky. Both maps correspond to a time-integration of 640 yr. 
The horizontal and vertical axes are given in arcseconds (for a 
distance of 5.8 kpc to the object). The Ha emission is depicted 
with a logarithmic scale (in units of erg s"' cm"^ sr"') given by 
the vertical bar on the right of the bottom panel. 



but for a jet propagating into the AGB wind (with a density 
distribution given by Eq.(|3). 

The gray-scale representations of the density stratifications 
obtained for integration times from 100 to 600 years (see Fig|2l 
andO show the head of the jet traveling away from the injec- 
tion point. The ejection velocity time variability produces inter- 
nal working surfaces that travel down the jet (which we iden- 
tify with the emitting knots of the observed jet). Each working 
surface is produced when fast material catches up with slower 
material ejected at earlier times. 

The temporal variability law used for the velocity in these 
simulations corresponds to the superposition of two terms : a 
periodic variability (described by a sinusoidal function) and 
a linear increase in the velocity. While a linearly increasing 
ejection velocity vs. time produces a single, accelerating lead- 
ing working surface (see Raga et al. 1990 and Canto, Raga & 
D'Alessio 2000), the sinusoidal mode (which is superimposed 
on the linear ejection velocity vs. time dependence, see Eq. 1) 
produces a train of internal working surfaces which travel down 
the jet beam (see, e. g., Raga & Noriega-Crespo 1998). 

A comparison between Fig|2land Fig |3l illustrates the ef- 
fect of the stratified environment on the evolution of the knots. 
Initially (at times < 300 yr), the head of the jet propagates 
faster in the uniform environment model. However, as the time- 
evolution proceeds, the head of the jet in the AGB wind case 
accelerates as it gets out to the lower density regions farther 
away from the source, and the jet eventually becomes longer 
than the one of the homogeneous environment jet model (see 
the t = 600 yr model of Figs.|2and|3. 

An interesting difference between the homogeneous envi- 
ronment and the AGB wind models can be seen in the Ha emis- 
sion maps (Figl^. While the uniform environment jet model 
shows extended bow shock wings (associated with the leading 
working surface), the AGB wind model does not show such 
wings. This is a direct result of the fact that at distances of 
~ 10"^ cm from the source the AGB wind has a very low den- 
sity, so that the emission from the bow shock is much lower 
than in the homogeneous environment model. 

In this way, we see that the model of a jet travelling into 
an AGB wind produces Ha intensity maps that more closely 
resemble the observations of the Hen 3-1475 jet, as this ob- 
ject does not show emission from extended bow shock wings 
associated with the jet head. Because of this, we carry out 3D 
simulations only for the case of a jet moving into an AGB wind. 



Fig. 5. Density (left), pressure (centre) and temperature (right) 
stratifications for a precessing jet moving into an AGB wind, 
obtained from the 3D numerical simulation described in the 
text. Each panel shows a cut (on the y - z plane) of the cor- 
responding variable for a 640 yr integration time. The scales 
of both axes are in units of 10'^ cm. The density, pressure and 
temperature scales are logarithmic (top horizontal bar on each 
panel), given in g cm""* dyn cm"^ and K, respectively. 

Fig. 6. [N II] 6583 A, [S II] (6717-H6730) A, [O I] 6300 
A,[0 III] 5007 A and Ha emission maps (from left to right), 
corresponding to an integration time of 640 yr. These maps 
have been obtained integrating the emission line coefficient 
along of lines of sight, considering that the jet axis has an angle 
of 50° respect to the plane of the sky. The horizontal and ver- 
tical axis are in arcseconds (after considering a distance of 5.8 
kpc) while the vertical bar on the right, shows (in logarithmic 
scale) the normalized emission. 

Fig. 7. HST image in [N II] of Southeastern outflow of Hen 
3-1475, which exhibits a jet-like structure and three pairs of 
symmetric emission knots. The innermost pair of knots show 
subarcsecond structure, with the existence of three well defined 
compact subcondensations. The knots have been labeled fol- 
lowing Riera et al. (2003). 

and do not study 3D models of jets moving into a uniform en- 
vironment. 

4.2. Comparison between 3D numerical results and 
HST observations 

Figure|5lshows cuts on the y-z plane of the density (left), pres- 
sure (centre) and temperature (right) stratifications, for the 3D 
simulation of a precessing jet moving into an AGB wind (see 
§3.2). These stratifications correspond to an integration time of 
640 yr 

From the density stratification (left panel of Fig|5j, we see 
that the jet has a leading head, followed by two high den- 
sity knots, which coincide with high-pressure regions (central 
panel). There is a high-temperature (> 10^ K) region which 
corresponds to the "cocoon", while the internal jet regions have 
temperatures ~ 10^ K or lower (right hand panel of Fig|5}- 

The main effect of the precession is to produce a curved 
morphology for the jet beam, which qualitatively resembles 
the "S"-shaped morphology of the Hen 3-1475 jet (see Fig|5}- 
Otherwise, the 3D model has internal working surfaces which 
are similar to ones of the axisymmetric models (see §4.1). 

The predictions obtained from the 3D model can be com- 
pared wiffi the WFPC2 HST images and with the position- 
velocity diagrams obtained from the (spatial and spectral) high- 
resolution STIS spectroscopy. In order to compare the 3D nu- 
merical simulation with the HST images of Hen 3-1475, from 
the model we have computed intensity maps for the Ha, [N II] 
6583 A, [O I] 6300 A, [S II] (6717-H6731) A and [O III] 5007 
A (from left to right. Fig. |6} integrating the emission coeffi- 
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cient along lines of sight. We have assumed that the y-axis lies 
in the plane of the sky, and that the z-axis is oriented toward the 
observer at an angle of 50° with respect to the plane of the sky. 
Furthermore, we have changed the physical coordinates of the 
model to angular sizes (after taking into account distance and 
projection effects). 

In Figure|6l we show the simulated line emission maps for a 
640 yr integration-time, which can be compared with the HST 
[N II] image of the SE jet (shown in Fig. 0. The images of 
Hen 3-1475 obtained with the HST using narrow band filters 
(such as F658N) essentially trace the ionized gas through the 
detection of the nebular emission coming from the correspond- 
ing emission line (such as [N II] 6583 A). However, we should 
note that losses are expected in the detection of the extremely 
blueshifted and redshifted emission due to the large velocities 
involved. Therefore, the observed intensities of the knots in the 
[N II] image are underestimated due to the reduced transmis- 
sion of the filter F658N. The losses are larger for the innermost 
knots, which show the largest radial velocities. 

All of the maps (from high to low-ionisation species) show 
the same number of knots (at the same distances from the in- 
jection point) with similar morphologies. Comparing the ob- 
served [N II] image (Fig. with the simulated emission line 
maps (Fig|6j, it is clear that the observed and simulated maps 
are qualitatively similar. 

Att = 640 yr (see Fig.|6j, we see the leading head at an an- 
gular distance of ~ 7'.'3 from the injection point, which agrees 
well with the observed angular distance of 7'.'57 from the cen- 
tral source to the outermost knot SE3. In all of the simulated 
emission line maps we observe that the leading head has a dou- 
ble structure. Actually, this double feature corresponds to two 
working surfaces. This pair of knots is formed by a faster work- 
ing surface which at f = 640 yr has just overtaken a slower 
working surface ejected earlier (this dynamical effect can be 
seen in the PV diagrams of Fig. 8). 

The intermediate knot (labeled SE2 in Fig.0 is located at 
a distance of ~ 5'.'9 in both the observed and simulated im- 
ages. The observed innermost knot (labeled SEla, SElb and 
SElc) is an elongated structure formed by three subconden- 
sations (at distances from the source from ~ 2'.'27 to 3'.'65). 
No emission is observed in the innermost region (i.e. between 
SEla and the central source). In the inner region, the simu- 
lated maps (see Fig. |6j show an elongated feature connecting 
the point of injection with a compact knot at ~ 4'.'5 down the 
jet. This feature clearly resembles the morphology observed in 
the innermost knot (SEl) of the Hen 3-1475 jet. However, the 
calculated emission maps do not show the observed emission 
gap between this knot and the source. 

A detailed comparison of the intensity maps of low- 
ionization (i.e. [O I] or [S II]) and high-ionization species (i.e. 
[O III]), reveals that the leading bow-shock is prominent in [O 
III] and Ha, and much fainter than the other knots in the low- 
ionization species (see, e.g., the [O I] intensity map of Fig|6}. 
The simulated [O I] intensity map shows a large range of inten- 
sities from the bright innermost knot to the fainter leading head. 
All knots show a more or less constant [O III] intensity (from 
the inner knots to the head). Unfortunately, no spectrophoto- 
metric observations of the three knots of Hen 3-1475 (cov- 



ering the wavelength range from 5000 to 6800 A) have been 
published. 

Figure|8ldisplays position-velocity (PV) diagrams in [Nil] 
6583 A (left) and Ha (right), obtained from our simulations for 
an integration time of 640 yr. Both maps have been convolved 
with a spectral resolution of 50 kms in order to simulate 
the STIS spectral resolution. In order to have predictions that 
can be directly compared with the spectroscopic observations 
(Borkowski & Harrington 2001; Riera et al. 2003; Riera 2004), 
the line profiles have been computed assuming a broad spectro- 
graph slit (covering all of the computational domain along the 
y-axis) straddling the z-axis. Again, an angle of 50° between 
the z-axis and the plane of the sky has been assumed. 

The calculated PV diagrams (see Fig. |8} qualitatively re- 
produce several characteristics of the observed kinematics. As 
expected, the large radial velocities and the observed decrease 
of the radial velocity with distance along the Hen 3-1475 jet 
are well reproduced by the model. Also, the models show dis- 
tinctive step-like radial velocity decreases along the jet axis, 
and such an effect is also seen in the observed PV diagrams 
(see Fig. Eland Riera 2004). 

The observed emission line profiles arising from the knots 
show wide double-peaked profiles (Borkowski & Harrington 
2001; Riera et al. 2003) . In the innermost knots the high- 
velocity peak is stronger than the low-velocity peak, while the 
situation reverses at the intermediate knots. The calculated PV 
diagrams show broad emission line profiles, but which are nar- 
rower than the observed profiles by a factor of ~ 2. Also, the 
calculated PV maps show double-peaked profiles only at very 
specific spatial locations (e. g., in the knot at x 4'.'7 from the 
source in Fig.|8}. 

5. Discussion and Conclusions 

We have carried out 2D (axisymmetric) and 3D numerical sim- 
ulations for modeling the jet of the proto-planetary nebula Hen 
3-1475, which is the first attempt for obtaining a full and dy- 
namical model for the outflows associated with this nebula. 

We have found that it is possible to fit the observed structure 
of the Hen 3-1475 jet with a time-dependent jet model with an 
ejection velocity history composed of a sinusoidal mode super- 
imposed on a linear ramp (i.e., a linear increase of the ejection 
velocity with time). The sinusoidal mode has a period of ~ 120 
years and a half-amplitude of 150 km s and the linear term 
has a slope of 1 km s"' yr"'. We have explored models of jets 
moving into a homogeneous environment, as well as models of 
jets moving into a previously ejected, stratified AGB wind. 

We should note that Lee & Sahai (2003) have presented 2D 
hydrodynamical simulations of the interaction of a CFW/jet 
and a spherical AGB wind in a PPN. Lee & Sahai adopted a 
periodic variation in the density and velocity of the jet (with 
an amplitude of 150 km s ' and a period of 22 yr, keeping 
a constant mass loss rate) in order to reproduce the morphol- 
ogy and velocities observed in the PPN CRL 618. A model 
including both periodic density and ejection velocity variations 
could describe the characteristics of PPNs such as Hen 3-1475. 
However, Lee & Sahai (2003) point out some difficulties for 
producing the knotty structure observed along the axis of PPN 
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Fig. 8. PV diagram for emission in [Nil] (left panel) and Ho- 
(right panel), for an integration time of 640 yr An angle of 
50° between the plane of the sky and the jet direction has been 
considered. The horizontal axis is given in units of km s ' while 
the vertical axis is in arcseconds. The [Nil] and Ha emission 
are shown in logarithmic scale by the bar on the right, in units 
of erg s"' cm"^ sr"' 

CRL 618. It is unclear whether or not this kind of velocity 
and density variability would work well for modelling the knot 
structure of the PPN Hen 3-1475. 

From the numerical models we have obtained simulated 
[Nil], [SII], [OI], [OIII] and Ha maps, as well as PV diagrams, 
which can be directly compared with the observational results 
which have been previously obtained with the HST. A compar- 
ison between the simulated emission maps (obtained from the 
2D and 3D numerical simulations) and the observations shows 
that a model of a variable speed precessing jet (with a period 
of 1500 yr and the ejection velocity history discussed above) 
moving into a previously ejected AGB wind produces the best 
agreement (of all of the computed models) with the observed 
characteristics of the Hen 3-1475 jet. 

The observed string of shock-excited knots is well repro- 
duced in our simulated maps. Also, the knot positions are sim- 
ilar to the observed ones, and are distributed along a curved 
locus which resembles the S-like shape observed in HST im- 
ages of the Hen 3-1475 jet (see Figs.|6landE}. 

With respect to the kinematical characteristics, by compar- 
ing the calculated and the observed PV diagrams we can say 
that our model succesfully explains : 

- the large observed radial velocities, 

- the decrease of the radial velocity of the knots with increas- 
ing distances from the source, 

- step-like changes in radial velocity at the positions of the 
knots, 

- the existence of broad emission line profiles. 

However, the simulated profiles from 3D simulations are 
narrower than the observed ones (by a factor of 2), and do not 
show double-peaked profiles with strong emission at the high- 
est velocities. 

From this discussion, we conclude that many of the mor- 
phological and kinematical features of the Hen 3-1475 jet 
can indeed be reproduced (at least in a qualitative way) by a 
model of a variable ejection velocity, precessing jet. This quite 
promising result gives a framework for modelling and inter- 
preting future, more detailed observations of this jet, which 
could lead to more stringent constraints for the proposed 
model. Finally, if taken seriously, the parameters of our model 
give an idea of what are the variability periods and amplitudes 
(as well as the precession) associated with the ejection from 
Hen 3-1475, which in principle provide constraints on the na- 
ture of the stellar source. 
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